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Abstract 

We consider the modified f{R) theory of gravity whose higher or- 
der curvature terms are interpreted as a gravitational fluid or dark 
source. The gravitational collapse of a spherically symmetric star, 
made up of locally anisotropic viscous fluid, is studied under the gen- 
eral influence of the curvature fluid. Dynamical equations and junc- 
tion conditions are modified in the context of f{R) dark energy and by 
taking into account the expansionfree evolution of the self-gravitating 
fluid. As a particular example, the Skripkin model is investigated 
which corresponds to isotropic pressure with constant energy density. 
The results are compared with corresponding results in General Rel- 
ativity. 
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1 Introduction 

Modifled f{R) theory of gravity constitutes an important development in 
modern cosmology and theoretical physics. One of the main motivations 
of this theory is that it may lead to some interesting results about dark 
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energy (DE). The other motivation comes from the fact that every unification 
of fundamental interaction exhibits effective actions containing higher order 
terms in the curvature invariants. This strategy was adopted in the study of 
quantum field theory in curved spacetimes jl] as well as in the Lagrangian 
of string and Kaluza-Klein theories [2] . 

Higher order terms always give an even number as an order of the field 
equations. For example, term produces fourth order field equations [3], 
term R\3R (where □ = V^V^) gives sixth order field equations P^, l5], simi- 
larly, R\3'^R yields eights order field equations [6] and so on. Using conformal 
transformation, the term with second derivative corresponds to a scalar field. 
For instance, fourth order gravitational theory corresponds to Einstein the- 
ory with one scalar field, sixth order gravity corresponds to Einstein gravity 
with two scalar fields, etc. [HIT]. In this context, it is easy to show that f{R) 
gravity is equivalent to scalar tensor theory as well as to General Relativity 
(GR) with an ideal fluid [8j. 

Let us now show that how f{R) gravity can be related to the problem of 
DE by a straightforward argument. When the Einstein-Hilbert (EH) gravi- 
tational action in GR, 

Seh = 77- / d'^x^/^R, (1.1) 
2k J 

is written in the modified form as follows 

SfiR) = d'xV^fiR), (1-2) 

the addition of a non-linear function of the Ricci scalar demonstrates to cause 
acceleration for a wide variety of f{R) function, e.g., [9]-[l8]. Variation of 
f{R) action with respect to the metric tensor leads to the following fourth 
order partial differential equations 

FiR)R^p~^fiR)g^l,-VaVpFiR)+g^^DFiR) = kT^P, («, /3 = 0, 1, 2, 3), 

(1.3) 

where F{R) = df{R)/dR. Writing this equation in the form of Einstein 
tensor, it follows that [H] 

Ga, = j{T-+T;f^^), (1.4) 
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where 



f{R) + RF{R) ^^^ ^ V, V;3F(/?) - g^pDF{R) 



:i.5) 



It is clear from Eq. (ll.4p that curvature stress-energy tensor T^^^ formally 
plays the role of a source in the field equations and its effect is the same as that 
of an effective fluid of purely geometrical origin. In fact, this scheme provides 
all the ingredients needed to tackle the dark side of the universe. Thus 
curvature fluid can play the role of both dark matter and DE. The "Dark 
source" term is not restricted to hold the usual energy conditions. Therefore, 
f{R) theory may be used to explain the effects of DE on cosmo logical and 
gravitational phenomena of the universe. 

The importance of gravitational collapse lies at the center of structure 
formation in the universe. A starting smooth arrangement of matter will 
eventually collapse and make the powerful structures such as stellar groups, 
planets and stars. One can study the gravitational collapse by taking the in- 
terior and exterior regions of spacetime. The proper junction conditions help 
to study the smooth matching of these regions. In Einstein gravity, the pio- 
neer work |20] was carried out on dust collapse by taking static Schwarzschild 
in the exterior and Friedmann like solution in the interior spacetime. 

During fluid evolution, self-gravitating objects may pass through phases 
of intense dynamical activities. The dynamical equations are used to observe 
the effects of dissipation over collapsing process. Skripkin [21] studied the 
central explosion of the spherically symmetric fluid distribution under con- 
dition of constant energy density. The evolution of such a fluid yields the 
formation of a Minkowskian cavity within the fluid distribution centered at 
the origin. This problem is studied in detail by Herrera et al. [221 123] • It was 
shown that under Skripkin conditions, the expansion scalar vanishes which 
requires the existence of a cavity within the fluid. Sharif et al. investigated 
some aspects of gravitational collapse regarding singularity and event hori- 
zons in f{R) theory [21] and also, in GR dissipative fluid collapse with and 
without adding charge [25]" EZ]- In a recent paper [2B], we have studied the 
effects of f{R) DE on a dissipative collapse. Some special solutions are also 
discussed for nondissipative case. 

In this work, we investigate that how DE generated by curvature fluid 
generally affects the dynamics of viscous self-gravitating fluid. In fact, this 
provides the generalization and extension of our previous work [2H]- We 
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explore the contribution of f{R) DE in developing Skripkin model. The 
plan of the paper is as follows. In next section 2, we formulate the field 
equations in metric f{R) gravity for spherically symmetric distribution of 
anisotropic fluid. Section 3 is devoted to derive junction conditions and 
dynamical equations with the inclusion of curvature fluid. In section 4, fluid 
evolution is investigated under the assumption of vanishing scalar expansion. 
Section 5 is devoted to study the Skripkin model in f{R) gravity. The last 
section 6 concludes the main results of the paper. 

2 Fluid Distribution and the Field Equations 

We consider a 3D hypersurface an external boundary of the star, which 
divides a AD spherically symmetric spacetime into two regions named as 
interior and exterior spacetimes. For the interior region to E^^\ we take the 
line element in the form as 

dst = A'it, r)de - B^t, r)dr'' - C^t, r){d0'' + sin^ 0d(f)''). (2.1) 

For the exterior spacetime to '£^'^\ we take the Schwarzschild spacetime given 
by the line element 

dsl ^ (1 - —)du^ + 2drdv - r\de^ + sin^ Odcj)''), (2.2) 
r 

where m represents the total mass and u is the retarded time. 

In the interior region, we assume a distribution of anisotropic self-gravitating 
fluid which undergoes dissipation in the form of shear viscosity. The energy- 
momentum tensor for such a fluid is defined as 

Tap = (P + P±)UaUi3 - P±gal3 + {Pr " P±)XaXl3 " '2Wal3- (2.3) 

Here, we have p as the energy density, p± the tangential pressure, pr the 
radial pressure, 77 the coefficient of shear viscosity, Ua the four- velocity of the 
fiuid and Xa the unit four- vector along the radial direction. These quantities 
satisfy the relations 

u'^Ua^l, = X'^Ua^O (2.4) 

which are obtained from the following definitions in co-moving coordinates 

= A-^S^, x" = B-^^l. (2.5) 
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The shear tensor aab is defined by 



(2.6) 



where the acceleration and the expansion scalar are given by 

tta = Ua;l3U^, Q = M° • (2.7) 

The bulk viscosity does not appear as it can be absorbed in the form of 
radial and tangential pressures of the self-gravitating fluid. From Eqs.f l2.5( ) 
and (12. 6p . we have the following non- vanishing components of the shear tensor 



2 1 

(Til = ~ 3^^^' ^22 = 3^*^^' ^33 = 0-22 slu^ 9, 

where a is the shear scalar and is given by 

1 f B C 

From Eqs. (l2.5p and (12. 7p . we have 



(2.8) 



(2.9) 



«i = — r 



A' 
A'' 



a°'ar 



A' 
AB 



e 



1 B C 



(2,10) 



where dot and prime represent derivative with respect to t and r respectively. 
In Einstein frame, the field equations (ll.4p for the interior metric become 
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(2.13) 
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Here the components of dark fluid are obtained from Eq. Ol.Sp as follows 
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The Ricci scalar curvature is given by 
A" 



R 



B AB 

+ 
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3 Junction Conditions and the Dynamical Equa- 
tions 

Here, we develop equations that govern the dynamics of dissipative spher- 
ically symmetric collapsing process. For this purpose, we use Misner and 
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Sharp formalism [22]. The mass function is defined by 



m r) = -{l + 9^-C,C,) = 2 h + - 5^ 



(3.17) 



From the continuity of the first and second differential forms, the matching 
of the nonadiabatic sphere to the Schwarzschild spacetime on the boundary 
surface, T.^^\ yields the following results 



M{t,r) = m 



(3.18) 



and 
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The important reason for generating the junction conditions at the stellar 
surface is to study the dissipative evolution of the star. Using the field 
equations fl2.12p and f l2.13p in Eq.f l3.19p . we obtain 



Pr - -W 



AB 



52 



(3.20) 



In the next section, we would discuss dynamics with physically meaningful 
assumption of vanishing expansion scalar which describes the rate of change 
of small volumes of the fluid. The expansionfree fluid evolution should imply 
the formation of a vacuum cavity within spherically symmetric fluid distri- 
bution ^22j. Taking E^*^ [i stands for internal) to be the boundary surface 
of that vacuum cavity and matching this surface with Minkowski spacetime, 
we have 



0, -Pr - -W 



S(') -til 



+ 



01 



52 ■ AB 

The proper time and radial derivatives are given by 
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l_d_ 



(3.21) 
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where C is the areal radius of a spherical surface inside the boundary E*^^\ 
as measured from its area. The velocity of the collapsing fluid is defined by 
the proper time derivative of C, i.e., 



U = DtC = ^ 
A 



(3.23) 



which is always negative. Using this expression, Eq. (l3.17p implies that 



E = — = [l + U^+ 



/2 



(3.24) 



When we make use of Eqs.lEH]), ( l2J0ll and (Km in Eq.(l2l2]), we obtain 



E 
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01 
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The rate of change of mass in Eq. fl3.17p with respect to proper time, with 
the use of Eqs. fl2.11l) - fl2.14l) . is given by 



DtM 
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rp{D)\ rpiD)' 

^11 \u-E^^' 



AB 



(3.26) 



This represents variation of total energy inside a collapsing surface of radius 
C. In the case of collapse ?7 < 0, the terms inside the first round brackets 
increases the energy density through the rate of work being done by the 
effective radial pressure pr + The presence of dark fluid component 

shows the contribution of DE having large negative pressure. These terms 
appear with positive sign representing negative effect, hence decrease the rate 
of change of mass with respect to time. Now, it depends upon the strength of 
DE terms that they may balance the positive effect of effective radial pressure 
or overcome on them. Likewise, we have 
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00 
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(3.27) 



This equation describes how energy density and curvature terms influence 
the mass between neighboring surfaces of radius C in the fluid distribution. 
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Here the rate would decrease in the consecutive surfaces by the repulsive 
effect of DE. Integration of Eq.f l3.27p with respect to "C" leads to 



M = 47r / — 
2F 



P + 



rp{D) J J rp{D) 

-^00 _ 

E AB 



dC. 



(3.28) 



The dynamical equations can be obtained from the contracted Bianchi 
identities. Consider the following two equations 
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Using Eqs.dMD, (EDD, (E22D and (EMD, it follows that 



C 



,(3.30) 



(3.31) 



~[DtP + ^(3p + Pr + 2px)e + ^(pr - P± - 2r7or)a], (3.32) 

E 



T"^;/3Xa = ^I^c(Pr + -r/o-) + (p + p, + -r/a)a + 2(p,-p^ + 2r7(T)-. 

(3.33) 

The acceleration DtU of the collapsing matter inside the hypersurface is 
obtained by using Eqs.daHS]), (jS^S) and ^7I^ as follows 



DtU 



M An 
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Substituting a from Eq. fl3.34p into fl3.33p . it follows that 



{p + Pr + -r]a)DTU 



-{P + Pr + -W) 



M Atc 4 T] 



11 



c 



4 2 



(3.35) 



This equation shows the role of different forces on the collapsing process. The 
factor within the brackets on the left hand side stands for "effective" inertial 
mass and the remaining term is acceleration. The first term on the right 
hand side represents gravitational force. The term within the first square 
brackets shows how shear viscosity and DE affect the passive gravitational 
mass. The first two terms in the second square brackets are gradient of the 
effective pressure and effect of local anisotropy of pressure with negative sign 
which increases the rate of collapse. 



4 Expansionfree Evolution of Self-gravitating 
Fluid 

In this section, we take = and discuss expansionfree evolution of the 
self-gravitating fluid. In such a case, Eq. fl2.10p yields 

B--^c ^ ^^-^^ 

where gi{r) is an arbitrary function. Implication of vanishing B requires that 
innermost shell of the fluid should be away from the center of the collapsing 
sphere. This situation initiates the formation of a vacuum cavity at the 
center [22]. The physical meaning of this condition is explained by taking 
two different definitions of radial velocity of the fluid. Substituting Eq.(4.1) 
and value of T^jf ^ from Eq. flZT^ in f ETT^ . we get 
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Integration yields 



A = r^e" 

Tit) 



J { F AF >c . 



(4.3) 



Thus the interior metric becomes 



nit) 
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dr 
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dr"^ - C\t, r)ide^ + sin^ ed(f)'^). 



(4.4) 



The above metric represents spherically symmetric anisotropic fluid which is 
going under shearing expansionfree evolution. 

To proceed further, it is necessary to assume the condition of constant 
scalar curvature (i? = Rc), according to which F(i?c) = constant. In /(-R) 
theory, the case of constant scalar curvature exhibits behavior just like solu- 
tions with cosmological constant in GR. This is one of the reason why the 
DE issue can be addressed by using this theory. In this case metric reduce 
to the following 



ds' 



FJCC^ 



dt^ - ^dr^ - C\t, r)ide^ + sin^ 



(4.5) 



Here we have chosen (71 (r) = 1 without loss of generality. For this metric, 
the scalar curvature (I2.16P becomes 
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One can find such values of C and ti for which the condition of constant 
scalar curvature is satisfied. For constant scalar curvature, the components 
of dark fiuid f l2.15p reduce to 

A^VL ~{Dc) « „(D^ „(n„\ —C^VL 
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where 

The corresponding junction conditions (IS.lSp . ( 13.20P and ( ]3.2ip will take the 
form 

M{t, r) ^= m pr ^= 1^, 

OTT 

M(t,r)^=0, P. = ^. (4.9) 

OTT 

Using Eq.diSD in the field equations (l2TT]) -( l234l) . it follows that 

-2C'C" - 5C'C" + - 3-^ = j,{8np + Q), (4.10) 

Idcct + ^ = 0, (4.11) 



= -^(Sirp, - n), (4.12) 

FFCC'S \ C nj yc CC C C2 

= h8np^-Q). (4.13) 
Similarly, Eqs. fl3.26p - fl3.27p become 

implying that 

M = ^ r (p+—]c'dC. (4.16) 
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Further, Eq.f l4?T4|) yields 



^vrp, = ^^ + a (4.17) 
C^C 
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For Schwarzschild mass m, Eq. fl4.17p is fully consistent with the junction 
conditions Eq.f l4.9p . 

On the similar conditions, Bianchi identities are reduced to the following 

P; + (P + VrVj: + 2(p + 2p, - - -i^ = 0. (4,19) 

We would like to mention here that in GR for isotropic fluid, i.e., = P±, 
Eq.f l4.18]) implies that energy density p depends only on r. However, in f{R) 
theory it remains the function of both time and radial coordinate. Integration 
of Eq. fl4.10p after some manipulation yields 

1 T2-2m 1 n / 1 1\ 

Here T2{t) is an arbitrary integration function. Using Eqs. fl4.5p and f l4.20p 
into Eq. fl3.17p . we obtain 

ft f I l\ , , 

^2 = -T 7- - • (4.21) 



k \67r F, 



5 Skripkin Model 

In the Skripkin model, it is assumed that fluid has isotropic pressure [pr = 
p± = p) and constant energy density, i.e., p = po. Here the fluid is as- 
sumed to be at rest initially and after that there is a sudden explosion at 
the center keeping the Skripkin condition. It is noted here that in GR, the 
expansion scalar automatically vanishes for the Skripkin conditions by virtue 
of Eq. (13.32 p . However, in f{R) gravity, it will not vanish and we take expan- 
sionfree evolution in order to study Skripkin model. Consequently, Eq. fl4.20p 
becomes 

where 

k = (5.23) 
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Under Skripkin conditions, Eq. fl3.35p takes the form 



(po + p)DtU = -(po + p) 
Using Eq.(I523) in Eq.(^J2i>, we obtain 



'f7 



3k + 2 



_r2_ 



For isotropic pressure, Eg. (14. 9 p becomes 



P 



in 



Substituting Eg. (15. 26 p in Eg. (15. 25 p and integrating, we get 



T2 



k + 



3F, 



(5.24) 



(5.25) 



(5.26) 



(5.27) 



where ci is an integration constant. Using Egs. (l4.5p and (I5.22p into Eg. (l3.17p . 
the mass function becomes 



(5.28) 



The total mass of the configuration m is obtained by measuring M on the 
boundary E^^^ as follows 



n Cl 

M = — Cwe) - — = m, 
3Fr ' 2 



(5.29) 



Substituting value of ci from Eg. (15.291) into Egs. (15.271) and (I5.28p . we have 



r2 



M 



k + 



n 



k + f]{C'-C^,.,)+m, 



(5.30) 



(5.31) 



Applying the matching conditions (I3.2ip on the boundary surface of the 
cavity inside the fluid. Eg. (15. 31 p yields 



(5.32) 
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Using this equation in Eq.f l5.3Up . we obtain 



r2= + (5.33) 



Differentiating Eq. (l5.32p with respect to "t", we have 
implying that 



^s(=)^2('=) ~ C'2(i)C's(') (5.34) 



C 

A-^(,)=A^[i) and ?7s(e) = — f^C/sW- (5-35) 

This result can also be obtained by definition of \J on the boundary surfaces. 
Applying Skripkin condition of constant energy density, Eq. (l4.16p yields 

Differentiating with respect to "t", we have 

Substitution of Eq.f l5.37|) into Eq.f l4.14]) with isotropic pressure yields 

p = --^ + F~ ( 1 ~ TT ) = constant. (5.38) 
Comparing Eq. (15.381) with Eq. fl5.26p . we obtain 

po = ^. (5.39) 

This shows that Skripkin model is eliminated by the junction conditions. It 
is mentioned here that in GR, po = while in our case, it turns out to be a 
non-zero constant. 

Since pressure will be maximum on some spherical surface C = Cs (say), 
hence pressure gradient must vanish on this surface. Thus it follows from 
Eq. flTO]) that 

(p + p)(CV)-- A^^O. (5.40) 

The solution of this equation yields the surface Cs which divides the fluid 
into two regions with a positive and negative pressure gradients respectively 
in the inner and outer sides of the surface. 
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6 Summary and Conclusions 



The most important feature of f{R) gravity is concerned with the issue of 
DE. In this paper, we have studied the effects of f{R) DE on the dynamics of 
a self-gravitating spherically symmetric star. The star is made up of viscous 
anisotropic fluid distribution and dissipating energy in the form of shearing 
viscosity. This work extends and generalizes our recent paper [28]. As a 
special case, the Skripkin model is studied which is based on the simplest 
conditions of isotropic pressure and constant energy density. 

This study is devoted to the expansionfree evolution of the fluid collapse 
resulting in the formation of a vacuum cavity within the fluid distribution. 
We have developed junction conditions on two hypersurfaces. One is the S*^^\ 
separating the fluid distribution from the Schwarzschild spacetime while the 
other is E*^*) , the boundary of internal cavity within which we have Minkowski 
spacetime. 

It is found that DE arising from the curvature fluid affects the whole 
dynamics of the gravitational collapse due to its repulsive effect. For example, 
it decreases the rate of change of mass of collapsing sphere with respect to 
time and adjacent surfaces. The comparison with the corresponding GR 
yields more general results and hence provides extra degree of freedom. For 
example, in Skripkin model, the comparison with [22], [23] is given as follows 

• In GR, implication of Skripkin conditions on the dynamical equations 
yields expansionfree evolution, i.e., G = 0. However, no such relation 
is obtained in f{R) gravity. 

• In GR, isotropic pressure with nondissipation leaves the energy density 
as only a function of radial coordinate. Here, Eq. fl4.18p does not imply 
that energy density is time independent due to presence of curvature 
fluid (even in the simplest case of constant scalar curvature). 

• Matching of results on the boundary surface S*^^^ gives a non-zero con- 
stant value (based on constant scalar curvature) to Skripkin energy 
density po while it becomes zero in GR. 

It it worth mentioning here that Skripkin model is completely consistent 
with the junction conditions. Moreover, the expansionfree models with the 
appearance of a vacuum cavity may be used to study the formation of voids 
at cosmological scales [30j for any type of fluid. 
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